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The temperature dependence of the magnetic susceptibility, x(T), is investigated for 
one- dimensional interacting electron systems at quarter-filling within the Kadanoff- Wilson 
renormalization-group method. The forward scattering on the same branch (the <?4-process) 
is examined together with the backward (<?i) and forward (g^) scattering amplitudes on oppo- 
site branches. In connection with lattice models, we show that x(T) is strongly enhanced by the 
nearest-neighbor interaction, an enhancement that surpasses one of the next-nearest-neighbor 
interaction. A connection between our predictions for x(T) and experimental results for x(T) 
in quasi-one-dimensional organic conductors is presented. 

KEYWORDS: magnetic susceptibility, renormalization group, long-range Coulomb interaction, extended- 
Hubbard model, quarter-filling, one-dimension, (TMTTF) 2 X, (TMTSF) 2 X 



In low dimensional conductors, the influence of long- 
range Coulomb interactions, or the off-site interac- 
tion V's are known to lead to various interesting or- 
dered states, such as charge ordering 1 and spin-density 
wave (SDW) state coexisting with charge-density wave 
(CDW). 2 ' 3 At the present stage, however, the effect of 
V's is not clearly understood compared with that of the 
on-site repulsion U. The purpose of this letter is to ana- 
lyze the effect of V's on the quasi-one-dimensional con- 
ductor in the normal state by calculating the tempera- 
ture dependence of the magnetic susceptibility x(T) for 
the one-dimensional system at quarter-filling. The results 
shed some light on the magnetic properties of both nor- 
mal and ordered states in these kinds of materials. 

A lot of studies of xCO were devoted to the one- 
dimensional (ID) Hubbard model, i.e., for V's = 0. The 
Bethe ansatz gives us an exact solution, but x only at 
T = is available for the ID Hubbard model. 4 The tem- 
perature dependence of x can be extracted from numer- 
ical simulations but only at high temperatures. 5 For U 
smaller than the bandwidth, the renormalization-group 
(RG) approach gives results that quantitatively agree 
with both the exact solution at T = and the numer- 
ical solutions at high temperatures. 5 In the presence of 
V's, however, there is no exact solution, and the size of 
numerical simulations becomes exceedingly large. These 
techniques are then of limited use to investigate the effect 
of V's. 

In this letter, we use the RG technique to calculate 
x{T) for the full temperature range from zero to the 
bandwidth Eq in the presence of V's (> 0). Consider- 
ing the possible interactions for branches of right and 
left going electrons in the continuum limit, we take into 
account the forward scattering on the same branch (SB), 
denoted as the (^-process, with the backward (gi) and 
forward (172) scattering processes on opposite branches. 
This is the first derivation of the RG flows of the sets 
of coupling constants and x(^) at the one-loop level in- 



cluding the non-logarithmic channels for particles on the 
SB, which become important at finite-temperatures. 

Our main results are: 1) x(T) is noticeably enhanced at 
finite temperature by the nearest-neighbor repulsion (Vl) 
compared with the next-nearest-neighbor repulsion (V2); 
2) x(T = 0) is not simply proportional to the inverse of 
the ordinary spin- velocity v a ; 3) xCO i s reduced by a 
moderately large V2 due to the CDW fluctuations. 

We consider the quarter-filled extended-Hubbard 
model in ID with the Hamiltonian H = Ho + Hj, where 
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Here t denotes the intrachain hopping integral, c\" a as an 
annihilation (a creation) operator on the i-th site with 
spin a (=t, |). The density operators are = c\ a Ci^ 
and rii = In weak-coupling regions, Eq. (1) can 

be safely expressed as 



p,k : cr 



(3) 



where e p (k) — vp(pk — kp) is the linearized ID dispersion 
and vp being the Fermi velocity. The operator c^ k a an- 
nihilates (creates) an electron of spin a close to the Fermi 
point of the right k = +kp (p = +1) and left k = — kp 
(p = —1) branches. We focus on the scattering processes 
between electrons near the Fermi points and express the 
Hamiltonian in terms of the charge- and spin-couplings 
in the following form: 

Hi =- Y [9pPp(<l)P-p(-<l) + 9aS p (q) ■ S- p (-q) 

+ 94pPp(q)p P (-q) + g4aSp(q) ■ S*(-qj] , (4) 
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(a) Cooper 



(b) Peierls 




(c) pp-SB 



(d) ph-SB 





(e) pair vertex part 



Fig. 1. Diagrams at one-loop level for (a) Cooper and (b) Peierls 
channels, and (c) pp and (d) ph channels on the same branch 
(SB); (e) pcrturbative expansion of the pair vertex part z a [see 
eq. (17)]. The solid (dashed) line indicates electrons for right 
(left) going state, and the shaded square denotes respective cou- 
pling constants. 



where 



Pp(l) — 2 ^ c \.k+a.o C v.k.cn 
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are the charge and spin-density of the branch p, respec- 
tively. Here a is the vector of the Pauli matrices. The 
<74 CT -term is equivalent to (g4a/3)S p (q) ■ S p (q) due to the 
spin-rotational SU(2) symmetry of the original extended- 
Hubbard Hamiltonian. The coupling constants at 1 /4 fill- 
ing (kp — 7t/4:,vf — v2t) are given by 



g p = U + AV 1 +6V 2 , 
g a = -U + 2V 2 , 

g 4p = U + 4V 1 +4V 2 , 

g ia =-U. 



(7) 
(8) 
(9) 
(10) 



The lattice constant is taken as unity. 

Following the Kadanoff- Wilson RG technique, 7 each 
RG step consists of the partial integrations of the fermion 
degrees of freedom in the outer band-momentum shell, 
E (£)/2 ^ ep(k) > E {£ + d£)/2 for electrons and 
-E {£)/2 ^ e p {k) < -E {£+d£)/2 for holes. Here, E (£) 
is the renormalized bandwidth Eo(£) = E$e~ l with the 
initial bandwidth Eq — 2vpkp ~ 2t. At the one-loop 
level (Fig.l), RG equations for coupling constants at fi- 
nite temperature are given by 

^G p (£) = -\g p {£) [G ip (£) ~ G 4tr (£)} I SB (£), (11) 



-G a {£) = Gl{£)I CF {£) 



-G a {£) [G 4p (£)-G ia (£)]I SB (£), 



(12) 



d_ 

If 
d 



G ip {£) 
GM 



[G 2 p (£) - G Ap (£)G 4a (£)] I SB (£), (13) 
[3G 2 a (£) - G Ap {£)G Aa {£)] I SB (£), (14) 



where G„(£)'s are dimensionless couplings with the ini- 
tial conditions G„(0) = g v ji:vp. The quantities /cp(^) — 
tanh[_Eo(^)/4T] are the Cooper/Peierls outer shell con- 
tractions [Fig. 1 (a), (b)], which are cut off by the tem- 
perature, namely Icp(^) — 1 for Eq(£)/4 > T and 
Jcp(^) — otherwise. 

The function Is B (£) is the coefficient of the particle- 
particle (pp) and particle-hole (ph) contractions on the 
SB [Fig. 1 (c) and (d), respectively] given by 



IsB(£)d£ 



n rp Shel l 
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(15) 



where & p (k,uj n ) = [iuj n — e p (fc)] 1 is the non-interacting 
Green's function with a Matsubara frequency uj n . The 
SB channels become finite due to thermal particle-hole 
excitations, i.e., E (£)/A ~ T. 

By setting I SB (£) = and I CP (£) = 1, Eqs. (11)- 
(14) reproduce the previous RG equations which treat 
only the Cooper and Peierls contractions. 6 ' 7 Equations 
(11)-(14) are the full RG equations at the one-loop level. 
Here, "full" means they include not only the logarith- 
mic (Cooper and Peierls) channels, but also the non- 
logarithmic (SB) channel. By considering the contribu- 
tion of 7sb, the renormalization of g^ P Aa and their in- 
fluence on the charge- and magnetic-susceptibility is ob- 
tained. In Eqs. (11)-(14), they apparently imply that the 
spin- and charge-part are coupled together. However, the 
contribution of the SB channel, /sb(^), is small compared 
to that of the logarithmic term, Jcp {£) > so that the cou- 
pling of spin- and charge-part is very weak. Further, at 
the ground state, the present flow equations exhibit the 
spin-charge decoupling since Isb{£) = at T = 0. 

The uniform magnetic susceptibility x(^) can be cal- 
culated by introducing a Zeeman coupling between a 
source field h p (q) and the spin-density variable S p (q). 
By adding Hh = J2 q P z ah p (q)S p (q) to the Hamiltonian, 
the uniform magnetic susceptibility can be obtained and 
which leads, 

2 f°° 

X(T) = / [z a {£)f I SB (£)d£, (16) 

where z a {£) is the pair vertex part and is calculated from 
-G a {l) + G^{i)]I SB {i), (17) 



with 



G 4± (£)= [G ip {£)-G ia {£)], 



(18) 



and the initial condition z a (0) = 1 [see Fig. 1 (e)]. 

The flows of the couplings as a function of energy at 
T = O.Oli are shown in Fig. 2. We show results for the 
following different sets of parameters; (i) (U, Vx,V 2 ) = 
(2i,0,0), (ii) (U,V U V 2 ) = (2M,0) and (iii) (U,V U V 2 ) = 
(2t, 0, i). Note that for the cases (i) and (ii), the SDW 
becomes dominant, while both the SDW and CDW are 
enhanced for the case (iii). Generally, for Eq(£)/A > T, 
G a {l) follows G a (£) = Go- (0)/ [1 + G„(0)ln(E (£)/E )}, 
as found in the conventional RG. With decreasing Eq(£), 
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Fig. 2. Energy flow of the couplings at T = O.Oli. (The SB 
channel becomes effective around Eq(1) ~ AT = 0.0l8i?o.) Go- 
and G4x(= (G^p — Gi a )/2) are shown by lines (i) and (i') for 
(U,V 1 ,V 2 ) = (2t,0,0), by lines (ii) and (ii') for (U,V 1 ,V 2 ) = 
(2t,t,0) and by lines (iii) and (iii') for (U, Vi, V 2 ) = (2t,0,t). 




Fig. 3. Magnetic susceptibility x(T) for (i) (U, Vi, V 2 )=(2t, 0, 0), 
(ii) (2i,i,0), (iii) (2t,0,t) and (iv) (2i,0, 2i), where xo = 
2/irvp(~ 0.45/t) denotes the magnetic susceptibility without in- 
teractions at T = 0. 



|G CT |'s [plots (i) & (ii)] are gradually suppressed due to 
the irrelevance of G a (< 0) in the limit of T = 0. How- 
ever, the G4j_'s [plots (i'), (ii') & (iii')] keep their initial 
values until Eq(£) decreases to Eo(£)/4 ~ T. Around 
Eq(£)/4 ~ T, the SB channel becomes effective while the 
contributions of the Cooper and Peierls channels to the 
flow are strongly suppressed. As a result, the SB chan- 
nel slightly enhances \G a \ and remarkably reduces Gi± 
[plots (ii) & (ii')] for finite V x . For finite V 2 , G 4 ± [plot 
(hi')] is reduced further. The contribution of the SB chan- 
nel is relatively small for the Hubbard model (Vi,2 = 0) 
[plots (i) & (i')]. As shown later, the properties of x(T) 
are directly affected by the temperature dependence of 
G a and G±±. 

The temperature dependence of the magnetic suscep- 
tibility x(T) is shown in Fig. 3. The result for the Hub- 
bard model (U — 2t) given by the plot (i) is quanti- 
tatively consistent with x{T) calculated by the quan- 
tum Monte Carlo (QMC) simulations. 5 The value at 



T = 0, x(0)£ = 0.58, agrees well with the exact solu- 
tion, Xcxact(0)t = 0.57, for U — 2t. 4 Near zero tem- 
perature, x(T) decreases as dx/dT — > oo for plots (i) 
and (ii), which has also been pointed out in the con- 
text of the Heisenberg spin chain. 9 This logarithmic de- 
crease of \(T ~ 0) is due to the logarithmic decrease of 
|Go-|. 5 Note that for large U, e.g., U > it, there was 
a sizable difference between the previous x(T) of the 
RG approach 5 and the x(^) °f the QMC calculation at 
high temperatures. The difference is much reduced by 
the present RG calculation due to taking account of the 
non- logarithmic terms (the SB channels), which become 
important at high temperatures. It is expected that such 
a difference is further reduced by the modified-Kanamori 
approach, 5 which is not applicable at low temperatures 
but gives better agreement with xqmc(T) at high tem- 
peratures. In the present calculation, we take the lin- 
earized (continuum) dispersion, in which lattice effects, 
coming from the use of the full tight binding spectrum 
c p (k) = —It [cos(fc) — cos(kp)), are neglected. When the 
lattice effects are considered in the calculation of Isb (£) , 
there appears a slight dip around T ~ O.K. 5 For U = IVi 
[plot (iii)] corresponding to the boundary between SDW 
and CDW, G a is invariant due to G a (0) = 0, so that 
dx/dT is constant at low temperatures as for an ordi- 
nary paramagnetic metal. For U < 2V2 [plot (iv)], x is 
reduced to zero at T <C t due to the CDW state. 

It is noticed that x(T) is much enhanced by Vi com- 
pared with V2 at finite temperatures. This difference orig- 
inates from the degree of the renormalization of . As 
seen in Fig. 2, in the presence of V\ , the reduction of 
G4± [plot (ii')] due to the SB channel is smaller than 
that in the presence of V2 [plot (iii')]. Consequently, the 
enhancement of x(T) by V\ is much larger than by V2 at 
finite temperatures. This superiority of V\ is seen only 
at finite temperature. Such an effect of V\ is surprising 
since V\ was considered to be irrelevant whereas V2 to be 
relevant. 8 

The origin of this enhancement can be further analyzed 
by considering the following random phase approxima- 
tion (RPA) for X (T) 



X(T) = 



xgen 



irvp 1 + [G CT (T) - G 4i _(T)]x£(T)/2' 



(19) 



where x2(T) = tanh(S /4T) is the magnetic suscepti- 
bility of the non-interacting case per branch. In Eq.(19), 
G a u±_)(T) denotes the renormalized coupling at T, i.e., 
G a (T) = G„(0)/[l + G a (0)ln{2T/E )] and G 4 ±(T) = 
[G 4p (£ = 0) - G 4a (e = 0)]/2. Note that this formula in- 
cludes the fluctuation beyond the simple RPA due to the 
renormalization of the Cooper and Peierls channel above 
T. The result of Eq. (19) qualitatively agrees with that 
of Eq. (16) except that it does not consider transients 
due to the interplay between the Cooper (or Peierls) and 
SB channel for E (£) - 4T. At T = 0, Eq. (19) yields 
x(0) = 2(itv*)~ 1 where v* is given by 

v* =v F [l- G 4 ±/2] 

= v f [1-{ttv f )- 1 {U/2 + V 1 +V 2 )]. (20) 
As is clear from this form, x(T) depends not only on U 
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Fig. 4. Contour plot of x(T)/xo at T = O.OOlt as a function of V\ 
and Vi. The plot at Vi/t = U/2t = 1 corresponds to the phase 
boundary between SDW and CDW. 

but also on V\ and Va- In earlier studies, x(0) is given 
merely by the inverse of v a , where the spin velocity v a is 
expressed as v a = Vp[l — U/(2itvf)] and is independent 
of V\ and V 2 . 6 The present results reveal that x(0) is not 
given by the ordinary spin- velocity v a , but by v*. Such a 
difference originates from a careful treatment of H\ [Eq. 
(4)] (especially the 54-process) and the accurate deriva- 
tion of the flow of z a [Eq. (17)]. It is worth remarking 
that in the rotationally invariant Hartree-Fock theory, 
the V\ interaction does contribute to the magnetic sus- 
ceptibility. The corrections to the spin velocity obtained 
in this work are consistent with those of the Hartree-Fock 
theory. We note that a recent result obtained by the nu- 
merical diagonalization also indicates the enhancement 
of x (0) by Vi. 10 

Figure 4 shows the contour plot of x m the V1-V2 plane 
at T = O.OOlt and U = 2t. It exhibits x( T ) is mono- 
tonically enhanced by increasing Vi.. For large V2, x(^) 
decreases exponentially to zero with increasing V 2 indi- 
cating the formation of CDW state. 6 ' 8 The V1-V2 depen- 
dence of x(T) is approximately proportional to V\ + Va 
for small V^/t. The magnitude of V 2 /t(~ 1.3) for vanish- 
ing of x is larger than that of T = 0, i.e., V% = U /2(= t), 
due to the effect of finite temperature. 

The role of V\ and V2 on x can be intuitively un- 
derstood as follows. First, consider the Hubbard model 
(U > 0, Vi = V2 = 0), where the dominant state is the 
2fc F -SDW with a periodicity of four lattice spacing. The 
enhancement of x(^) by U(> 0) occurs due to the re- 
duction of double occupancy of up spin and down spin. 
On the other hand, at low temperatures x(T) is reduced 
due to irrelevance of the backward scattering. 6 Next, the 
addition of V\ , which prevents electrons from coming to 
the nearest- neighbor site, enhances the 2fc F -SDW and 
the local moment of spin. This explains the enhance- 
ment of x(T) in the presence of V\ as seen from plots (1) 
and (2) in Fig. 3. When V2 is added, electrons become 
less located at next-nearest-neighbor sites leading to the 
coexistence of 2fc F -CDW and 2fc F -SDW. The CDW fluc- 
tuations reduce xCO- 

Finally, the present results can be discussed in con- 
nection with the experimental situation for the series of 



quasi-one-dimensional organic conductors (TMTCF) 2 X 
(C=S, Se; X=PF 6 , Br). 11 In (TMTTF) 2 X for example, 
x(T) exhibits a maximum and dx/dT is large at low tem- 
peratures. The present results suggest that besides the 
effect of U, a moderate V\ leads to a clear maximum 
in x{T) and to large dx/dT. Considering that a charge 
ordered state exists in (TMTTF^X, 1 these compounds 
should be characterized by a large V\ . Therefore the pre- 
vious determinations of U from X {T) in (TMTTF) 2 X 
might be considered as overestimations. 11, 12 As for the 
Bechgaard salts, no charge ordered state is found in 
(TMTSF)2PF6 and a much less pronounced maximum of 
x(T) and smaller dx/dT are observed. Besides the role of 
U, the present results would indicate that (TMTSF) 2 PF 6 
has a relatively large V 2 (and not small U and V\). A 
large V2 is consistent with a 2fcp-SDW state coexisting 
with 2fc F -CDW. 2 - 3 

In conclusion, we calculated the magnetic susceptibil- 
ity, x{T), by considering the renormalization of the for- 
ward scattering in the same branch, the f/4-process, us- 
ing the Kadanoff- Wilson RG technique. With the care- 
ful treatment of the (74-process, we found that x(^) is 
enhanced by the nearest-neighbor interaction V\, but is 
less enhanced by V2 • It is also found that x(0) is given by 
v* [Eq. (20)], which is different from the ordinary spin- 
velocity v a . A comparison with the experimental situ- 
ation of the quasi-one-dimensional organic compounds 
suggests that (TMTTF) 2 X has a large Vi (< U) but 
small V 2 , while (TMTSF) 2 A has a large V 2 (< V X <U). 
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